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ABSTRACT 

o : 
o . 

CN ' Ultraviolet absorption from interstellar '^CO and ^^CO was detected toward p 

^ • Oph A and x Oph. The measurements were obtained at medium resolution with 

^ . the Goddard High Resolution Spectrograph on the Hubble Space Telescope. Column 

^ : density ratios, N{^'^CO)/N{^'-^CO), of 125 ± 23 and 117 ± 35 were derived for the 

I sight lines toward p Oph A and x Oph, respectively. A value of 1100 ± 600 for the 

^ ; ratio A^('^C^^O)/iV(^^C'^0) toward p Oph A was also obtained. Absorption from 

' vibrationally excited H2 {v" = 3) was clearly seen toward this star as well. 

■ The ratios are larger than the isotopic ratios for carbon and oxygen appropriate for 
Tjj- ■ ambient interstellar material. Since for both carbon and oxygen the more abundant 

. isotopomer is enhanced, selective isotopic photodissociation plays the key role in the 

. fractionation process for these directions. The enhancement arises because the more 

I abundant isotopomer has lines that are more optically thick, resulting in more self 

Qh| shielding from dissociating radiation. A simple argument involving the amount of self 

O I shielding [from A^(^^CO)] and the strength of the ultraviolet radiation field premeating 

' the gas (from the amount of vibrationally excited H2) shows that selective isotopic 

■ photodissociation controls the fractionation seen in these two sight lines, as well as the 
^ ■ sight line to ( Oph. 

H 
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Subject headings: ISM: abundances - ISM: molecules - stars: individual (p Oph A, x 
Oph) 

1. Introduction 

Carbon monoxide is the second most abundant molecule, after H2, in interstellar clouds 
and is seen in dark cloud envelopes through observations at ultraviolet wavelengths (Smith & 
Stecher 1969; Morton 1975). High-quality spectra acquired with the Goddard High Resolution 
Spectrograph (GHRS) on the Hubble Space Telescope allow detailed studies of the CO abundance 
and the relative abundances of the various forms containing carbon and oxygen isotopes (Sheffer 
et al. 1992; Lambert et al. 1994; Lyu, Smith, &: Bruhweiler 1994). For instance, Lambert et 
al. deduced the relative abundances of ^^C^^O, ^^C^^O, ^^C^^O, and One goal of these 

studies is to extract information on physical conditions for the gas probed by CO absorption. The 
earlier efforts with GHRS of Sheffer et al. (1992), Lambert et al. (1994), and Lyu et al. (1994) 
had as a focus the sight line toward C Oph [RA(2000) = 16'* 37"* 10^; DEC(2000) = -10° 34' 02"]; 
here, we present results on two other sight lines, p Oph A [RA(2000) = 16'* 25"* 35*; DEC(2000) 
= -23° 26' 49"] and x Oph [RA(2000) = 16'* 27"* 01*; DEC(2000) = -18° 27' 22"], in the same 
portion of the sky. 

Prom observations of dark cloud cores via CO emission (e.g., Penzias 1981; Langer & Penzias 
1993) and of diffuse clouds via CH"*" absorption (Stahl et al. 1989; Centurion Sz Vladilo 1991; 
Crane, Hegyi, & Lambert 1991; Stahl &: Wilson 1992; Vladilo, Centurion, & Cassola 1993; 
Centurion, Cassola, &: Vladilo 1995), ambient isotopic ratios for C and O can be obtained 
because CO in cloud cores and CH+ in diffuse clouds are not fractionated. [Throughout this 
paper we consider (chemical) fractionation as any process that alters isotopic ratios from ambient 
values.] In cloud cores the severe attentuation of ultraviolet radiation removes all possible routes 
for fractionation, while the non-thermal conditions leading to CH"*" production are believed 
to equilibrate CH^" isotopomcrs. The ratios found from these studies of material in the solar 
neighborhood are: ^^C/^^C ~ 65, ^^O/^^O ~ 500, and ^^O/^'^O ~ 2600. These ^^C/^^C and 
16q^18q pa^jQg aj.g consistent with the recommendations of Wilson & Rood (1994) who claimed 
H2CO gives a lower limit, and CO an upper limit. Wilson & Rood obtained respective average 
values of 77 ± 7 and 560 it 25. The measurements reported here also probe material near the Sun. 

In cloud envelopes, however, the isotopic ratios in CO are altered by chemical fractionation. 
Isotopic exchange reactions (■'^^C"'" -|- ^^CO ^^C"*" -|- ^^CO) can enhance the abundance of ^^CO 
relative to -"^^CO when gas temperatures arc low (Watson, Anicich, & Huntress 1976) because 
-^^CO has a lower zero-point energy. On the other hand, the more abundant variants are enhanced 
relative to less abundant ones when selective isotopic photodissociation dominates (Bally & Langer 
1982; Chu & Watson 1983). Since CO photodissociation takes place via line absorption, the more 
abundant forms of CO have lines that are more optically thick, resulting in self shielding against 
further photodissociation for these isotopic variants, van Dishoeck and Black (1988), Kopp et al. 
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(1996), and Warin, Benayoun, &; Viala (1996) incorporated these processes into detailed models, 
and here we compare the predictions of these models with our observations. 

In our earlier studies of CO fractionation in the diffuse clouds toward ( Oph (Sheffer et al. 
1992; Lambert et al. 1994), significant enhancement of ^^C^^O relative to ^^C^^O was found, with 
a ratio about a factor of 2 larger than the interstellar value of 65 - 70. Lyu ct al. (1994), from 
the same data set analyzed by Sheffer et al., derived a much smaller ^^C^^O-to-^^C^^O ratio that 
was indistinguishable from the ambient value. One potential cause for the different conclusions 
comes from the set of oscillator strengths (/-values) for the A — X system of bands used to extract 
column densities. The analysis of Lambert et al. (1994) was based on the measurements of Chan, 
Cooper, & Brion (1993) from electron-energy-loss spectra, while the work of Lyu et al. (1994) 
adopted the /-values of Eidelsberg et al. (1992) from absorption measurements with a synchrotron 
source. Lambert et al. (1994) noted that their data with a higher signal-to-noise ratio yielded a 
more satisfactory curve of growth when the results of Chan et al. (1993) were utilized. Recent 
laboratory results (Smith et al. 1994; Federman et al. 1997a; Jolly et al. 1997, Zhong et al. 1997; 
Stark et al. 1998; Eidelsberg et al. 1999), based on a variety of experimental techniques including 
absorption of synchrotron radiation, are consistent with those of Chan et al., and therefore the 
appropriate set of /-values to use is no longer an issue. Furthermore, Lambert et al. derived 
i2(~ii6Qy'i2Qi8Q g^j^jj i2(-ii6Qyi2Qi7Q p^tios that Were also about a factor of two greater than 
the ambient ratios for the oxygen isotopes, a result expected if their conclusions about carbon 
fractionation were correct. Still, it is unsettling that two groups reached such different conclusions 
about CO fractionation, and one of the reasons for acquiring the observations toward p Oph A 
and X Oph was to compare the predictions of van Dishocck & Black (1988) for diffuse clouds with 
differing physical conditions. As described below, this objective was met. 

Ratios of the CO isotopic variants in diffuse and translucent clouds have been obtained via 
millimeter-wave techniques as well. Langer, Glassgold, &: Wilson (1987) mapped a region around 
Oph in ^^CO and presented a tentative detection of ^^CO toward the star. Their ^^CO/^'^CO ratio 
of about SOtlo is formally consistent with our earlier ultraviolet measure. More recently, Kopp et 
al. (1996) mapped the gas in the vicinity of ^ Oph in both isotopomers and obtained lower limits 
of about 45 to 65 for the more diffuse directions. Gredcl, van Dishoeck, Sz Black (1994) surveyed 
millimeter-wave emission from southern translucent clouds with Ay between 1 and 4 mag. While 
the ^^CO/"'^^CO ratios of antenna temperature, which may be affected by optical depth in the 
^^CO line, are low, the ^^C^^O/^^C^^O and ^^C^^O/^^C^^O ratios are similar to or greater than 
the ambient interstellar values. The ratios involving the oxygen isotopes suggest that selective 
isotopic photodissociation is operating in these southern clouds. 

It is not a simple matter to compare results from ultraviolet absorption and millimeter-wave 
emission lines. Absorption samples an infinitesmal pencil beam while emission is an average 
over the finite telescope beam. Furthermore, emission lines are more prone to uncertainties in 
excitation, radiative transfer, and abundance. Liszt & Lucas (1998) overcame these difficulties by 
measuring millimetric absorption against background compact extragalactic sources. They found 
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a ■'^^CO/-'^^CO ratio between 15 and 54 that decreases with increasing iV(^^CO). Such a result 
indicates that isotope exchange prevails in their sample of relatively diffuse clouds with iV(^^CO) 
as large as 2 x 10^^ cm~^. 

The paper is organized in the following manner. The next Section provides details of the 
ultraviolet measurements acquired with GHRS as well as ground-based observations of CH+, which 
provided the ambient ^^C/^^C ratio. Section 3 presents our results while Section 4 describes our 
analysis. In the latter section, results on the excitation of the fine structure levels in the ground 
state of C I and on absorption from vibrationally excited H2 are included as an aid in constraining 
the appropriate physical conditions for the clouds. The data on C I and H2 were extracted from 
the spectra used for CO absorption. The resulting column densities for C I appear in Zsargo, 
Federman, & Cardelli (1997), who derived a self-consistent set of C I /-values for future analyses. 
[The recent update by Federman &; Zsargo (2001) does not alter these column densities because 
column densities and Doppler parameters were inferred from lines with well-known /-values.] The 
relative populations in the fine structure levels are used to derive estimates for gas density and 
temperature (e.g., Jenkins, Jura, & Loewenstein 1983; Lambert et al. 1994) in the present work. 
Estimates of the flux of ultraviolet radiation permeating the gas are possible from analyses of 
vibrationally excited H2 (Federman et al. 1995; Meyer ct al. 2001). Next, simple arguments 
highlight the importance of selective isotopic photodissociation for the three Sco-Oph sight lines. 
General comparisons with the theoretical predictions of van Dishoeck & Black (1988), Kopp et al. 
(1996), and Warin et al. (1996) are made in the Discussion, where we suggest areas for further 
improvement in the theoretical models. Such improvements in our models (van Dishoeck &; Black 
1988) are beyond the scope of this paper. Throughout this paper we note specific isotopes as 
needed; otherwise, the general chemical formula is used. Similarly, explicit notation is given for all 
CO bands that are not part of the A — X Fourth Positive System, which is designated only by 
v'—O. An Appendix provides observational results for lines of elements heavier than zinc that are 
seen in our spectra. 

2. Observations and Data Reduction 

2.1. Ultraviolet Measurements 

Only a brief description of our reduction procedure is given here; details can be found in 
our earlier papers (e.g., Lambert et al. 1995). Spectra were acquired with grating G160M of the 
GHRS and covered three wavelength intervals - 1252 to 1294 A, 1385 to 1426 A, and 1443 to 1488 
A. Each exposure consisted of 4 FP-SPLITs and several exposures were obtained at each nominal 
wavelength setting, but slightly offset from each other. The most important step in the reduction 
involved the realization of the highest possible signal-to-noise (S/N) ratio from the data. We 
applied the procedures outlined by Cardelli &: Ebbets (1994) for the minimization of noise due to 
granularity, etc. on each FP-SPLIT subexposure. The final spectrum for a wavelength interval was 
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derived from merging the individual FP-SPLITs in wavelength space. Table 1 provides a listing 

of the exposures for each star, as well as other pertinent information, such as central wavelength 
and the number of counts. The total number of counts for an interval, the S/N ratio attained 
in the merged spectra, and the CO bands covered by the interval are also shown. The exposure 
times were 1632*, except for z2b9020g where one FP-SPLIT subexposure failed. The resultant 
S/N ratio in our final spectra is on average very similar to the tabulated values computed from 
(total counts)°-^. 

Two- Angstrom segments centered on each CO band were extracted and the stellar continuum 
was rectified with routines available in the IRAF environment distributed by NO AO. The 
rectification process still retained residuals that constitute the largest errors in the profiles of 
the strong bands {v' = 2-5). While uncertainties in equivalent width (Wx) as large as 25% are 
possible, the self-consistent column densities for the suite of bands suggests smaller systematic 
errors. A further problem arose while rectifying the spectra of x Oph. Since x Oph is a Be star, 
its continuum varies much more than do those of p Oph A and ( Oph. This introduced the largest 
source of error in the band profiles for the star because the acceptable rectified continuum depends 
to some extent on the person performing the task. Our final continua reflect the consensus of at 
least two rectifiers (Y.S. and S.R.F.). Figures 1 and 2 show respective spectra revealing absorption 
from ^^CO and ^^CO toward our targets. The observations are indicated by filled circles, and 
the fits (described below) to extract the ^^CO/^^CO ratio are shown as solid lines. Table 2 
displays the W\ values for the CO bands derived from the fully reduced spectra of p Oph A and % 
Oph. These include the intersystem bands, a' — X 14—0 and e — X 5—0. The root-mean-square 
(rms) variation in the continuum and the number of pixels across a band were employed in a 
conservative computation of the uncertainty in Wx. The CO bands in the spectra of p Oph A are 
much stronger (optically thicker) than those of x Oph and Oph. 

Absorption from C I and vibrationally excited H2 were seen in our spectra. The derivation 
of Wx for lines of neutral carbon toward both stars is presented in Zsargo et al. (1997). For 
absorption from v" = 3 in the lowest electronic state of H2, we employed a strategy much like the 
one discussed above for CO. Absorption involving the R(0), R(l), R(2), R(3), P(2), and P(3) lines 
was clearly detected in the spectrum of p Oph A (see Figure 3). The presence of the P(2) line 
manifested itself by an asymmetry in the line at 1279.478 A of C I from J = 2. Vibrationally 
excited H2 was not detected toward x Oph. The results for measures of Wx appear in Table 3, 
which includes the results of Federman et al. (1995) for 1^ Oph for reference. Our observations of 
all the expected lines in essence provide confirmation of vibrationally excited H2 in diffuse clouds, 
as suggested earlier by Federman et al. (1995). [The data of Meyer et al. (2001) are likely probing 
gas under more extreme conditions.] A steeply varying stellar continuum prevented the derivation 
of upper limits for R(0) and R(l) toward x Oph; furthermore, while the lack of an asymmetry in 
the C I line at 1279.478 A indicated that P(2) is not present, a useful upper limit could not be 
obtained. 
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2.2. Ground-based Measurements 

High-resolution observations of ^^CH+ and ^^CH"*" absorption toward p Oph A and x Oph 
were obtained at McDonald Observatory with an echelle spectrograph on the 2.7 m Harlan J. 
Smith Telescope in 1995 May. The stellar spectra were imaged onto a Texas Instruments CCD 
with 15 jim pixels. An interference filter was used to limit the spectral coverage to the single order 
containing the lines near 4230 A. Bias and flat-field images were acquired each night and Th-Ar 
comparison spectra were interspersed among the stellar spectra. Dark frames were taken the first 
night as a check on background levels during the longest exposure. In all, p Oph A {B = 5.26) 
was observed for 20 hours and x Oph {B = 4.70) for 14 hours. 

The raw data were reduced with standard IRAF procedures. After bias subtraction, flat 
fielding and removal of scattered light and cosmic rays, the spectra were extracted. Wavelength 
calibration was based on lines seen in the Th-Ar comparison spectra. The spectral resolution, 
which was derived from the full width at half maximum (FWHM) of the Th-Ar lines averaged 
over all the nights, was determined to be 33.1 it 0.3 mA, which is significantly larger than the 
thermal width for the Th lines (~ 2 mA). For p Oph A, fitting of the continuum level around the 
CH+ lines was straightforward. However, for x Opli) underlying stellar emission line caused 
some concern. Centurion et al. (1995; hereafter CCV) used a synthetic stellar line (which they 
assign to Fe II A4233.167) to extract and normalize the spectrum of the star HD110432 in the 
direction of the Southern Coalsack. An unfortunate fact of our spectrum was that only part of the 
stellar line was included, and so we performed a normal continuum division, restricted to a limited 
wavelength interval in order to minimize the influence of baseline structure. Since for x Oph the 
interstellar line is located on the blue wing of the stellar line, rather than in the self-reversal core 
(cf. CCV), we do not believe that our extraction procedure added significantly to the derived 
uncertainties. The rms dispersion in the rectified continuum yielded a signal-to-noise ratio per 
pixel of about 650 for both stars. The rectified spectra appear in Figure 4. 

After normalization/rectification, the W\ for the ^^CH"^ and ^^CH"*" lines were derived from 
Gaussian fits. Two methods utilizing the STSDAS routine NGAUSS were employed. In the 
first, we adopted an isotope shift of —0.272 A, as measured from the ^^CH'^ line. This value is 
based on the laboratory measurements by Bembenek (1997) for ^"^CH^ and those of Carrington 
& Ramsay (1982) for ^2CH+. We kept the widths of the ^^CR+ and ^'^CR+ hues the same in 
our syntheses. The results of this analysis appear in Table 4, where Wx, vlsr, and the Doppler 
parameter (6-value) are given. The LSR velocities obtained here are within 1 km s"^ of those 
derived by Crane, Lambert, &: Sheffcr (1995) from ultra-high-resolution observations; the 6-values 
also agree very well. The relatively large 6-values seen for CH+ lines are likely a consequence of its 
production involving an endothermic reaction (e.g., Lambert & Danks 1986). As for the W\, our 
determinations and those of Lambert & Danks (1986) are very similar, as are the results of Crane 
et al. for x Oph. For p Oph A, the measurements of Crane et al. yield an W\ about 2-a larger. 
The values of Wx yield column density ratios, N{^'^CR+)/N{^^CB.+), of 120 ± 54 and 65 ± 19, 
respectively, for the gas toward p Oph A and % Oph from a curve-of-growth analysis with / = 
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5.5 X 10~^ and b = 2.5 km s"-*^. Use of a smaller 6-value, 2 km s~^, introduces a negligible change 

because the optical depth at line center for ^^CH+ is less than 0.5. The uncertainties in the ratios 
are dominated by the precision of the ^^CH+ measures. The result for x Oph is in the range (60 
- 70) determined for other sight lines, but that for p Oph A is a bit high (consistent at the 1-a 
level, however). 

In the second method, we fixed the relative wavelength offset between the isotopomers at 
—0.265 A. This isotopic shift is based on the measurements of CCV for CH+ toward HD110432, 
rather than the value of —0.26 A obtained from theoretical calculations by Auguson and Herbig 
(1967) and used in earlier studies (e.g.. Centurion &; Vladilo 1991). The quoted uncertainty of 
about 2 to 3 mA in each laboratory measurement (Carrington Sz Ramsay 1982; Bembenek 1997) 
indicates that the shift derived from astronomical spectra is consistent with the laboratory value 
at about the l.5-a level. With a shift of -0.265 A we find 109 ± 36 and 61 ± 20 for p Oph A and x 
Oph. In summary, the ratios derived here indicate that the ^^C/^^C ratio is approximately 65, as 
in other, more precise determinations for the general interstellar medium (e.g., Langer & Penzias 
1993; CVV) and for the Rho Ophiuchi Molecular Cloud in particular (Bensch et al. 2001). 



3. Relative Abundances of CO Isotopomers 

Three steps were employed in the extraction of the relative abundances among the isotopic 

variants of CO. The extraction was based on synthesis of the measured bands through the 
minimization of the rms diff'erence between the observed band and the fit. The /-values of Chan 
et al. (1993) were adopted, as was the case in Lambert et al. (1994). The first step involved fits 
to absorption from the weaker 10-0 and 11-0 bands in order to set the initial column density. 
For this column density, simultaneous fits of the stronger bands yielded rotational excitation 
temperatures {Tj-ot = Tic T21, and T^2) and the 6- value. The final step involved fitting the ^'^CO 
bands, which are not very sensitive to changes in Trot and 6- value because these bands are optically 
thin. As consistency checks on these derived rotational excitation temperatures, Trot also was 
estimated from the 10-0 and 11-0 bands of ^^CO and the 2-0, 3-0, and 4-0 bands of ^^CO 
under the assumption that Tiq = T21 = r32. These determinations yield similar excitation 
temperatures for ^^CO, within the mutual uncertainties, and show that the rotational excitation 
in ^^CO and ^^CO is essentially the same. One added complication with the 2-0 and 3-0 bands 
for ^^CO was the presence of the C^^O band. For the synthesis of the ^^CO bands toward p Oph 
A, where all bands were much stronger than those toward x Oph, the /"^"^C^^O ratio was 

inferred as well. Figures 1 and 2 also show the synthetic spectra and the residuals between the 
observed and synthetic profiles. For completeness, the data are displayed with theoretical curves 
of growth, which are based on our results in Table 5, in Figure 5. All data can be described very 
well by these curves of growth, which were computed for a given N{CiO) with the appropriate 
wavelengths for the various lines comprising the bands. Allowance was also made for absorption 
from ^^C^^O in the strongest bands. 
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This procedure provided the necessary information for our study of chemical fractionation in 

CO; the results appear in Table 5. The values for N(^'^CO) are the most precise among available 
ones for the sight Uncs to p Oph A and x Oph. We note that the CO (A — X) bands toward p 
Oph A are stronger than the ones seen in the spectrum of ^ Oph (Sheffer et al. 1992; Lambert 
et al. 1994; Lyu et al. 1994). The ^^CO column density is not as large, however, because the 
6-values are greater. Our determination of the column density toward p Oph A is about a factor 
of 3 greater than the value inferred from Copernicus data of the B - X 0—0 band at 1150 A 
(Snow & Jenkins 1980). Most of the difference arises from the fact that they adopted a value of 
1.2 km s^"*^ where we deduced b = 0.6 km s^^. Because we fitted several bands simultaneously, 
whose wavelengths fall on different parts of the blaze function for G160M, our 6-values are less 
susceptible to the variations in spectral resolution across the blaze noted by Lyu et al. (1994). 
As for X Oph, the Copernicus measurements of Frisch (1980) on the C - X 0—0 band near 1088 
A, using a band oscillator strength of 0.123 (Fedcrnian et al. 2001), yield a value for A''(^^CO) 
that agrees nicely with ours. We also extracted the Copernicus data on p Oph A and x Oph for 
both the B — X and C — X 0—0 bands from the MAST archive at the Space Telescope Science 
Institute and used routines in IRAF for an improved comparison with the GHRS results. Table 
6 shows the comparison of measured and derived values of W\\ the derived values are based on 
the /-values of Federman et al. (2001) and the cloud parameters [iV(^^CO), Tio, T2\, and 6-value] 
from the synthesis of the A — X bands. The correspondence between the values is very good. The 
relatively slight differences for the C — X band are mainly the result of imprecise subtraction of 
the CI I line at 1088 A. The Wx for the C - X band toward p Oph A from Federman et al. (1980) 
is significantly less than that reported here; an incorrect continuum placement is likely the cause. 

While the above comparison shows that the derived cloud parameters are rather robust, we 
now comment on the velocities, 6-values, and excitation temperatures in Table 5. The heliocentric 
velocities associated with the weighted blend of the Q(l) and Q(2) lines for the 2-0 through 5-0 
bands agree with those obtained by Crane et al. (1995) to within about 2 km s~^. Since thermal 
velocities do not contribute significantly to the line width, we would expect similar 6-values for 
different species. Indeed, the 6-values of ~ 0.60 to 0.70 km s^^ agree with other measures based 
on ultra-high-resolution observations. Measurements of K I A 4044 toward x Oph (Knauth et al. 
2003) indicate a similar value. The K I column density from A 7699 (Welty & Hobbs 2001) can be 
made consistent with the results from A 4044 by assuming one velocity component with a 6-value 
of about 0.70 km s~^. Moreover, Welty & Hobbs suggest that the strongest component toward p 
Oph A has a 6-value of approximately 0.60 km s~^. For the direction toward p Oph A, synthesis 
of the strong and weak band profiles yielded different rotational excitation temperatures. While 
the 2—0 to 5-0 ^^CO bands indicate Tio of 2.7 K, the weaker 10-0 and 11-0 bands, as well as 
2-0 to 5-0 bands of ^^CO, are better fitted with Tio between 7.7 ± 2.5 and 10.3 ± 2.9 K. (The 
larger excitation temperatures are more consistent with T21 and T32 derived from the strong bands 
of ^^CO.) If the effect is substantiated by more precise determinations, it could be understood in 
terms of how deeply the observations probe into the cloud: optically thin bands, which sample 
most of the volume, reveal higher excitation temperatures in the core, but saturated bands are 
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sensitive to outside regions of the cloud. This effect can arise from a density gradient which also 
appears in our analysis of C I levels and vibrationally excited H2 described below. The effect of 
a decreasing CO excitation temperature toward the edge of the cloud can also result from less 
trapping of CO millimeter lines in the periphery compared with the amount of trapping in the 
cloud's core (e.g., Bernes 1979; Hogerheijde k. van der Tak 2000). An explicit radiation transfer 
code would be needed to confirm this scenario. In its stead, we also suggest that some combination 
of two clouds with differing physical conditions, such as the two dominant components seen by 
Welty &: Hobbs (2001) in K I absorption, might account for Tiq variations in the bands. Crane et 
al. (1995) found two CH components toward p Oph A as well, one with a small 6- value like that 
inferred by CO and the other with a larger 6-value. 

Two intersystem bands of CO, a' — X 14—0 and e — X 5—0, are present in our spectra of p 
Oph A and the stronger a' — X band is seen toward x Oph. Table 2 gives the measured values 
of W\. Federman et al (1994a) derived an /-value for the a' — X band which was (31 ± 7)% 
smaller than that computed by Morton k. Noreau (1994), while the /-value for the e — X band 
was (89 lb 16)% smaller. Fitting the bands with the parameters deduced from the A - X bands 
suggests f{a' - X) of (0.7 ± 0.2) /(MN) and /(e - X) of (1.1 ± 0.3) /(MN), where /(MN) is the 
value quoted by Morton & Noreau (1994). Our more precise determinations from spectra of X 
Per taken with the Space Telescope Imaging Spectrograph (Sheffer, Federman, k, Lambert 2002) 
reveal respective /-value ratios of 0.94 zt 0.15 and 0.78 =b 0.12. The combination of the results for 
these directions suggests that the band oscillator strength quoted by Morton k Noreau should be 
decreased some 10 - 30%. Improved analysis that includes multiple curve crossings (Rostas et al. 
2000; Eidelsberg k Rostas 2003) substantiates the astronomical findings. 

4. Analysis 

4.1. Density from Atomic and Molecular Excitation 

The relative populations of the fine structure levels in the ground state of C I yield estimates 
for gas density and temperature (e.g., Jenkins et al. 1983). The populations are mainly affected 
by collisional (de) excitation and far infrared radiative decay, and to a lesser extent by ultraviolet 
pumping and the subsequent decay from the excited electronic level. We previously analyzed the 
distribution of these levels in om study of the gas toward Q Oph (Lambert et al. 1994). Here 
we apply the same technique to the data on C I toward p Oph A and x Oph. Zsargo et al. 
(1997) derived column densities for each level; the resulting ratios, N{J')/N{J), are N{1)/N{0) 
= 0.430 ± 0.024 and 0.422 ± 0.019, respectively, for p Oph A and x Oph, and Ar(2)/iV(0) = 
0.238 lb 0.014 and 0.235 ± 0.011 for the two sight lines. The abundances (relative to total protons) 
of the collision partners (H, 0-H2, and P-H2) are also needed; these were obtained from Savage 
et al. (1977). For He, we assumed an abundance of 10%. As in all analyses of excitation, if the 
temperature is known, the density of collision partners {uc) is inferred. For neutral interstellar 
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gas, ric ~ n(H) + n(H2). 

The results of our analysis appear in Figure 6, where the allowed ranges in density and 
temperature are shown for 2-cr excursions about the observed column density ratios. One facet is 
immediately clear: The results for A'^(1)/A^(0) (dashed lines) suggest lower densities than do the 
results for A^(2)/A^(0) (dot-dashed lines). The earlier analysis by Jenkins ct al. (1983), based on 
spectra from the Copernicus satellite, did not reveal this dichotomy. Their column density ratios, 
though generally consistent with our estimates, probably had uncertainties too large to discern any 
difference. We also note that Jenkins et al. and Federman, Welty, & Cardelli (1997b) found this 
dichotomy toward the nearby line of sight to 0^ Sco, while Zsargo k. Federman (2003) detected the 
two velocity components in high-resolution ECH-A data on 1, 5, and a Sco. Ultra- high resolution 
spectra reveal additional, weaker components in K I absorption toward p and x Oph (Welty & 
Hobbs 2001). Moreover, for our two directions, the 6- values deduced from curves of growth for 
the lines originating from J = and 1 are larger than the 6-value for the J = 2 lines (Zsargo et 
al. 1997). From this information, we infer that absorption from the lower fine structure levels are 
probing (1) more extended portions of the neutral gas where the average density is expected to be 
lower or (2) an additional low density component. 

The use of the rotational excitation temperature for the J = and 1 levels in H2 [roi(H2)] as 
a measure of the kinetic temperature allows us to derive values for the gas density. Adopting the 
2-0" range in values of Savage et al. (1977), which arc also indicated in Fig. 6, yields densities for 
the gas toward both of our targets of « 100 and 200 - 400 cm'^ from N{l)/N{Q) and N{2)/N{Q), 
respectively. The estimates are rather robust. There is little change when applying substantial 
differences in the relative fractions of collision partners; the inferred densities are mainly influenced 
by the observed relative fine structure populations. For p Oph A, the effects of increasing the 
amount of UV pumping are illustrated in Fig. 6 as well. A 10-fold increase in UV flux, as 
suggested by our H2 measurements, lowers the density estimate a small, but noticeable amount. 

We can derive other estimates for gas density from analysis of the distribution of rotational 
levels in the ground state of C2 or CO. Excitation of C2 involves a combination of collisions, 
radiative decay between rotational levels, and near infrared pumping to the A electronic state, 
followed by radiative cascades among vibrational and rotational levels of the ground state. 
According to van Dishoeck & Black (1982), the amount of excitation can be represented by 
Uc ojlir^ where a is the coUisional cross section and Jj^ is the enhancement in IR flux over the 
typical interstellar value. The observational data on C2 of Danks & Lambert (1983) and van 
Dishoeck &: de Zeeuw (1984), along with the theoretical predictions of van Dishoeck & Black 
(1982), lead to the estimates shown in Figure 6 as enclosed areas. These were obtained through a 
X-squared minimization procedure, weighted by the precision of the observational data on column 
density for a given level. We adopted a cross section of 2 x 10~^^ cm~^, an oscillator strength of 
1 X 10~^, and Iir of 1. 

The results from C2 excitation are generally consistent with those from the N{1)/N{0) ratio 
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for C I, which is a bit of a surprise. The estimates from C2 are very sensitive to assumptions 
concerning excitation. Recent determinations (Langhoff et al. 1990; Erman & Iwamae 1995; 
Lambert et al. 1995) of the oscillator strength for the A — X 2—0 band, /20, are 20 to 40% larger 
than the adopted value. A larger /-value would increase the density estimate (van Dishoeck & 
Black 1982) a corresponding amount. The adopted cross section is in the middle of the range 
suggested by van Dishoeck & Black. Quantal calculations (Lavendy et al. 1991; Robbe et al. 
1992; Phillips 1994) span the suggested range, with larger cross sections giving lower densities. 
A reexamination of C2 excitation, utilizing the more recent determinations of /20 and collisional 
cross sections, appears necessary. 

Finally, the excitation of CO can be studied. This excitation is controlled by collisions, 
radiative decay, absorption of the cosmic background radiation, and resonant scattering of CO 
lines from nearby molecular clouds (Wannier, Penprase, &: Andersson 1997). The latter process 
appears to be the dominant one, and as a result, only upper limits on density are possible. This 
conclusion is based on the J = 1 — > maps of the Sco-Oph region by de Geus, Bronfman, & 
Thaddeus (1990). We selected the channel with vlsr between 3 and 5 km s~^ and estimated 
the filling factors for the gas at the projected positions of p Oph A and x Oph. The excitation 
temperature, Tiq, due solely to resonant scattering of emission from the molecular cloud is 7.3 ± 1.0 
and 3.2 it 0.4 K, respectively. These values are comparable to, or greater than, the values inferred 
from profile synthesis of the UV absorption bands. Therefore, all that can be said about density 
is that is it less than about 1000 cm~^. The fact that Tio is not larger than T21 or T32, as would 
be expected for subthermal conditions in relatively diffuse gas, is further proof that resonant 
scattering of line radiation dominates CO excitation along the two sight lines. 



4.2. UV Flux from H2 Excitation 

In cold (T < 1000 K) gas, the population of the v = 3 level of the ground electronic state is 
governed by the flux of ultraviolet radiation permeating the cloud (Black &: van Dishoeck 1987). 
Absorption from the J = 0, 1, and 2 levels in the vibrational ground state - the ones with the 
greatest populations - of ultraviolet photons leaves the H2 molecules in excited electronic states, 
but rapid decay returns the molecules to various vibrational levels of the ground electronic state 
(e.g.. Black & Dalgarno 1976). Collisional (dc)excitation is not effective in (de) populating v = 3 
because the densities and temperatTircs are too low. Thus, from the amount of absorption we can 
infer the ultraviolet flux (Federman et al. 1995). 

The measurements of vibrationally excited H2 toward p Oph A are not easily modeled. The 
amount of H I along the line of sight is very large, which, together with the large amount of H2 
in V = 3, suggests that there is low density gas close to the star. [Use of the revised H I column 
(4.3 X 10^^ cm~^) of Diplas Sz Savage (1994) in the analysis does not remove the difficulties 
encountered here.] On the other hand, the CO column density is quite high and Tio(H2) is only 
46 K (Savage et al. 1977), indicating the presence of denser, lower temperature gas along the 
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line of sight. Rotational excitation in C2 is consistent with a total density of a few hundred 
cm~^, although the quality of the data is not very high (see Fig. 6). [The total density, ntot = 
n(H I) + 2n(H2) > nc, is the parameter of interest in chemical modeling. For our sight lines, ntot 
is about 1.5 Uc-] 

The results from a set of models for p Oph A appear in Tabic 7, where the first five entries 
are the input parameters for the models: the enhancement factor for the UV flux (/(yy), n-jot, 
the gas temperature (T), the polytropic index, and the means of populating the levels during H2 
formation. All our models are based on an adopted 6-value of 1.0 km s~^. The initial model had 
constant ntot and T with ntot = 300 cm~^, T = 45 K, and luv = 4. This model reproduces the 
H2 = observations quite well, but it fails to reproduce those for H I and H2 f = 3. Increasing 
the radiation field to luv = 10 brings the column density for H2 in t; = 3, J = into agreement 
with observations, but this fails for the higher rotational levels. The H2 levels v = 3, J = 1-3 
clearly indicate a much higher rotational temperature, between 100 and 500 K; the second model 
shows the effects of increasing the temperature to 100 K. This improves the comparison with 
the rotational levels in v = 3, but now the column densities ior v = are not consistent with 
observations. An attempt based on polytropic models with a temperature gradient from about 
400 K at the edge to 45 K or 100 K in the center still could not reproduce all the data. Given 
the limited set of data, especially the lack of data for the high-lying rotational levels in v = 0, we 
can only state that there seems to be at least two components along the line of sight. A warm 
component of low density gas is exposed to a radiation field that is enhanced by a factor of ~ 10 
over the standard radiation field (Drainc 1978); presumably, this gas is located close to the B2 
V star, about 0.2 pc away. As noted in our work on ( Oph (Federman et al. 1995), the current 
models do not include line overlap, which can result in an overestimate of the pumping rate by up 
to a factor of 2. Thus, the luv inferred here for this component should be regarded as a lower 
limit. The second component, containing most of the H2 v = 0, as well as CO and C2, must be 
much colder and denser. This picture suggests that we are observing a photodissociation region 
in the reflection nebula associated with p Oph via absorption lines. The two components may 
correspond to the dominant components in K I (Welty & Hobbs 2001) and the two seen in CH 
spectra (Crane et al. 1995). Both K I and CH probe material over a larger range in density than 
either CO or C2. 

The direction to x Oph was modeled by van Dishoeck &: Black (1986); their best model (G) 
yielded predictions for the amount of H2 in f = 3 which exceed our upper limits by factors of 
2-4 or so (see Table 7). Therefore, the value for luv used in the calculations must be reduced 
accordingly, from which we infer that the strength of the ultraviolet field is comparable to the 
average interstellar field {luv of 1 to 2). As shown in Table 7, new calculations with reduced 
field strengths confirm this supposition. The results of the simplified analysis of Federman et al. 
(1994b) now agree better with the refined models described here. 
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4.3. CO Fractionation 

With the addition of the present results for p Oph A and x Oph, there are three sight Unes in 
See 0B2 with accurately determined values for the N{^'^CO)/N{^'^CO) ratio. We have 125 ± 23 
for p Oph A, 117 ± 35 for x Oph, and 167 ± 15 for ( Oph (the latter from Lambert et al. 1994). 
The amount of fractionation is comparable in the gas toward p Oph A and x Oph and it is larger 
toward C, Oph. These ratios are greater than the ambient ^^C/^^C ratio of about 65; self shielding 
of the more abundant isotoponicr clearly causes the observed fractionation. The results can be 
understood in a qualitative way through the realization that larger column densities increase the 
effects of self shielding, while larger fluxes of ultraviolet radiation increase the photodissociation 
rate. In other words, the 'reduced column density' NQ-'^CO)/Iuv is a measure of the amount of 
fractionation. Taking A^(^^CO) from Table 5 and Iijy from the previous section, along with the 
results for C Oph from Lambert et al. (1994), N{^'^CO) = 2.5 x 10^^ cm , and Federman et al. 
(1995), luv = 1 to 2, wc obtain values of N{'''^CO) / luv of about 1.9 x 10^^, 1.9 - 3.8 x 10^^, and 
1.3 — 2.5 x 10^'^ cm^^ for p Oph A, x Opli) and C, Oph, respectively. The 'reduced column density' 
follows the trend seen in the observed amount of fractionation. For the directions considered 
here, selective isotopic photodissociation is the process controlling the relative abundances of the 
^^CO and ^^CO isotopomers; the present results on C^^O are not precise enough for definitive 
statements. 

A more quantitative means of comparison involves analysis based on eqn. (15) in Lambert et 
al. (1994): 
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In this expression, the amount of fractionation, F13, is [nQ-'^CO) /nQ-^ CO)]/ [n{}'^C) /ni^'^C)] and 
the photodissociation rate for isotopomer z is Fj. The numerical coefficients in square brackets are 
based on the rate constants for isotopic charge exchange, while represents the elemental carbon 
abundance relative to the value for the Sun. A factor of 0.4 more closely matches the interstellar 
measurements of Cardelli et al. (1993a). The photodissociation rates come from the calculations 
of shielding functions by van Dishoeck &: Black (1988) — see their Table 5, which includes the 
effects of CO self shielding and mutual shielding from H2. The necessary columns of H2 are taken 
from Savage et al. (1977). The amount of grain attenuation is based on Model 2 described by van 
Dishoeck & Black (1986). For p Oph A and x Oph, we multiplied the grain optical depth by 0.7 
to accommodate larger than average values for the ratio of total to selective extinction (Federman 
et al. 1994). We approximate local densities n(X) with column densities. We considered Ijjv of 1 
toward x and C, Oph and luv = 10 toward p Oph A. For ntot-, we adopted 300 cm~^ for the gas 
toward p Oph A and % Oph, as inferred above, and 200 cm~^ for the gas toward C, Oph (Lambert 
et al. 1994). Finally, since van Dishoeck k, Black (1988) considered a slab geometry, the column 
densities used here are one half the measured ones. 
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With the use of an ambient isotope ratio of 65 for carbon, our measurements indicate values 
for Fi3 of 1.9 ± 0.4, 1.8 ± 0.6, and 2.6 ± 0.3 for the directions p Oph A, x Oph, and C Oph. The 
uncertainties in the observed values are based on the uncertainties in the ^^CO/^^CO ratios and 
an assumed 10% uncertainty in the ^^C/^^C ratio, taken in quadrature. The quantities in square 
brackets are approximately 1 when ntot of a few hundred cm~^ is adopted, as found from C I 
and C2 excitation. The main exceptions are the calculations for % and C, Oph where the bracket 
involving is about 2.0. In other words, this re-enforces the notion that the three sight lines 
are controlled by photodissociation. We find F12, to be 1.8, 1.0, and 1.2 for the respective sight 
lines. Considering the simplications made (e.g., column for local density), the accuracy of the 
observational results (10-30%), and the results using Table 5 of van Dishoeck & Black (1988) 
(20-30%), the agreement between calculations and observations is good. The correspondence for 
X and C, Oph could be improved by combined changes in luv and ntot of about a factor of 4; this 
would be accomplished by increasing luy and decreasing n^^. It is satisfying that our simple 
analyses utilizing observational input are able to reproduce the enhanced ^^CO/^^CO ratios seen 
by us. A more detailed modeling effort is beyond the scope of this paper. 

5. Discussion 
5.1. H2 Excitation 

When the populations of vibrationally excited H2 are controlled by optical pumping via 
absorption of UV photons, the ortho (odd J) to para (even J) ratio in excited states differs from 
the ratio for the vibrational ground state (Sternberg & Neufeld 1999). While their work focuses on 
warm photodissociation regions, where the ortho to para ratio in the ground state is the thermal 
value, the same ideas apply to cool gas where the population in J = is large. The effect is clearly 
seen in our data for p Oph A and the data of Federman et al. (1995) for C, Oph. The respective 
ortho to para ratios for the -y = 3 state toward the two stars are 2.7 =b 0.9 and 4.1 it 2.2, and the 
corresponding ratios for the ground state are 0.7 ± 0.3 and 1.2 ± 0.4. For the i; = 3 state, we 
summed over all possible rotational levels, including statistical weights, but only considered the 
J = and 1 levels of the ground state because they have much larger columns than higher-lying 
levels. 

The difference in ortho to para ratio between vibrational levels arises from enhanced optical 
depth in lines from J = of the ground state. Since optical pumping populates the rotational 
levels of vibrationally excited states via a two-step process, the change in rotational quantum 
number is or 2. The larger optical depth of lines originating from J = 0, the most populous 
(paxa) level, reduces the abundance of para rotational levels in v = 3. This leads to a larger ortho 
to para ratio for v = 3. 

Excitation of H2 also occurs during its formation on grains. During formation about 4.4 eV of 
binding energy is distributed among substrate excitation, molecular excitation, and translational 
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energy. The partition among the channels is poorly understood. For instance, it is not clear 

whether the nascent molecules are vibrationally hot and rotationally cold — high v and low J 
(see Duley & Wilhams 1993) - or the opposite (Hunter k Watson 1978; Wagenblast 1992). 
More recent experimental (Gough et al. 1996) and theoretical (Kim, Ree, &: Shin 1999) work on 
H2 excitation during formation on carbonaceous surfaces suggests most of the energy goes into 
vibrational motion. Internal excitation of H2 could be discerned from observations of high-lying 
rotational levels in the ground state and of vibrationally excited H2. Since high-lying rotational 
levels in the ground state are needed to assess the importance of the predicted distributions 
(Fcdcrman ct al. 1995), but such data do not exist for p Oph A, more definitive statements 
cannot be made at this time. The large number of ro-vibrational levels probed by observations of 
HD 37903 (Meyer et al. 2001) may provide the answer to this important question. 

In our studies of molecular chemistry in diffuse clouds (e.g., van Dishoeck &: Black 1989; 
Federman et al. 1997c; Knauth et al. 2001), we have begun to appreciate the importance of 
using the extinction curve for the sight line in the analysis. As for the two directions examined 
here, the extinction curve for both stars is lower than 'typical' at the shortest wavelengths 
(Fitzpatrick &: Massa 1990; Snow, Allen, & Polidan 1990; Green et al. 1992). In the present work, 
we did not attempt to model more than H2 in any detail. Under these circumstances, different 
extinction curves could be incorporated into the scaling term luv- It is also worth mentioning 
that calculations adopting bi-directional radiation fields, like those of Lyu et al. (1994), with one 
half the flux incident on each side of the cloud yield results similar to the ones presented here. 



5.2. CO Fractionation 

There are now three sight lines, all in Sco-Oph, where ^^CO is enhanced relative to other 
isotopomers by approximately a factor of two over the ambient ^^C/^^C ratio. We are confident 
in our results for a number of reasons. First, other analyses that yielded smaller relative amounts 
of ^^CO from the CO data toward (" Oph were based either on a set of /-values now known to be 
inaccurate (Lyu et al. 1994) or on data of the same bands for ^^CO and ^^CO (Levshakov & Kegel 
1994) which are prone to optical depth corrections in the more abundant isotopomer. The former 
point is addressed in more detail below. Second, our synthesis code does not always give results 
indicating enhanced amounts of ^^CO. In a paper on photodissociation regions toward sources of 
reflection nebulae (Knauth et al. 2001), our code used for profile synthesis confirmed the lower 
^^CO/^^CO ratio ('--' 40) seen toward 20 Aql (Hanson, Snow, & Black 1992) from an independent 
extraction of lU E spectra. 

An especially critical one is the third reason. New measurements of /-values for Rydberg 
transitions of CO are larger than once thought (Federman et al. 2001). [These are the /-values 
giving excellent fits to the CO bands seen in Copernicus spectra, upon adoption of the parameters 
derived from the A — X bands, as discussed in Section 3 and in Federman & Lambert (2002).] 
Although only a few transitions of importance in CO photodissociation have been re-examined to 
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date, in all cases the /-values for Rydberg transitions are about a factor of two larger than those 
compiled by Eidelsberg et al. (1991) from the measurements of Letzelter et al. (1987). This may 
explain why ab initio modeling efforts, such as those of van Dishoeck & Black (1986), van Dishoeck 
& Black (1988), Kopp et al. (1996), and Warin et al. (1996), are unable to produce ^^CO/^^CO 
ratios much in excess of the ambient ^^C/^^C ratio. Use of a larger /-value would increase the 
optical depth of dissociating transitions, which in turn leads to more self shielding for the more 
abundant isotopomer. Moreover, larger optical depths in lines leading to CO dissociation would 
improve the correspondence between model predictions and submillimeter-wave observations of 
dense photodissociation regions (HoUenbach & Ticlens 1999). Further work is needed to derive 
accurate /-values for all transitions involved in CO photodissociation. 

More self shielding for ^^CO would also provide a better correspondence in the CO abundance 
between model and observations. For instance, the theoretical predictions of van Dishoeck & 
Black (1986) are too small by factors of a few, even after accounting for the smaller values for Ijjv 
derived from vibrationally excited H2. An enhancement in CO production would help alleviate 
the problem as well. A key reaction is C+ -|- OH C0+ -|- H. Dubernet, Gargaud, & McCarroU 
(1992) performed quantal calculations on this reaction and obtained a rate constant of about 
5 X 10^^ cm'^ s~-^, which is substantially larger than the estimate given in Prasad &: Huntress 
(1980) and is about a factor of 2 larger than the value given by Federman &: Huntress (1989) based 
on average dipole orientation (ADO) theory. 

While a self-consistent picture seems to be emerging, we briefly discuss the apparent 
contradictory results of Lyu et al. (1994) and Liszt k. Lucas (1998). It is our contention that the 
differences between Lyu et al. and our earlier work (Sheffer et al. 1992; Lambert et al. 1994) for 
C Oph lie in the adopted /-values. As noted in the Introduction, Lambert et al. found a more 
satisfactory curve of growth with the /-values of Chan et al. (1993). All agree that the column 
density of ^^CO is best determined by the 11—0 and 12-0 bands, which are essentially optically 
thin. For these bands, the /-values of Eidelsberg et al. (1992) are 50% larger than those of Chan 
et al. With their preferred curve of growth, Lyu et al. infer NQ^'^CO) = 1.8 x 10^^ cm ; this 
value becomes 2.7 x 10^^ cm~^ when using Chan et al.'s /-values, much more in line with the 
results of Lambert et al. (2.5 x 10^^ cm~^). We also point out that the two sets of results for the 
5-0 band of ^^CO agree nicely. There is another way to check the ^^CO/^'^CO ratios: the ratio of 
/-values for ^^CO and ^^CO bands with the same W\. Lambert et al. provided a comparison in 
their Table 5, which indicates a ratio greater than 125. The same comparison can be made with 
our data displayed in Table 2. For both p Oph A and x Oph, the Wx for the 11—0 band of ^^CO 
is the same as that for the 4—0 band of ^^CO. The ratio of /-values is 134, consistent with our 
isotopomeric ratios. 

Finally, we comment on the work of Liszt & Lucas (1998), where the ^^CO/^^CO ratio 
decreases with increasing Ar(-'^^CO) for column densities sampled by our UV measurements. First, 
our results specifically concern the clouds in Sco-Oph; other clouds, like that toward 20 Aql, 
have much smaller ratios (Hanson et al. 1992; Knauth et al. 2001). Second, there are no direct 
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measures of extinction or C"*" abundance for their extragalactic sight hnes. Enhancements in ^^CO 
relative to ^^CO require a significant flux of UV radiation to produce the C"*" needed for isotopic 
charge exchange (but not too large to heat the gas via photoelectrons coming off grains). A more 
accurate determination of the extinction would also help clarify why complex molecules such as 
HCN (e.g., Lucas & Liszt 1994) are detected, when standard chemical models of diffuse molecular 
gas (e.g., van Dishoeck & Black 1986) suggest very low abundances. 

In summary, ultraviolet observations of CO and its isotopic variants toward p Oph A and x 
Oph were analyzed in conjunction with measurements on absorption from C I and vibrationally 
excited H2. Significant enhancements in the amount of ^^CO relative to the other isotopomers 
were found. The amount of H2 in t; = 3 toward p Oph A indicates an ultraviolet flux about 
10 times the average interstellar flux permeates this direction; for the x Oph sight line, the flux 
cannot be greater than twice the average interstellar value. Analysis of C I (and C2) excitation 
reveals modest densities for the material seen in absorption. Resonant scattering of emission lines 
from nearby molecular clouds controls the distribution of CO rotational levels, thereby limiting 
the usefulness of CO excitation for inferring density. Simple arguments based on the processes 
involved in selective isotopic photodissociation are able to reproduce the observed ^^CO/^^CO 
ratios. Incorporation of an improved set of /-values for dissociating transitions in large scale 
models will likely lead to better agreement with our observations. 
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A. Abundances of Heavy Elements 

In addition to lines of CO, H2, C I, S I, Co II, and Ni II reported here and in Zsargo et al. 
(1997), Mullman et al. (1998), and Zsargo &: Fcdcrman (1998), our spectra of p Oph A and x Oph 
reveal absorption from Ga II, As II, and Sn II. The results are given in Table 8, where W\ and 
A^(X) are listed. The sources of /-values used in deriving column densities are Morton (1991) for 
Ga II, CardeUi et al. (1993b) for As II, and Schectman et al. (2000) for Sn II. A curve of growth 
with a 6-value of 2.5 km s~^ was used; such a 6-value is common for the dominant ion (Savage, 
Cardelli, &: Sofia 1992). Two components, separated by about 10 km s~^, are discerned in the 
spectra of x Oph. These are the major complexes found in Na I D (Welty, Hobbs, & Kulkarni 
1994). Furthermore, knowledge of the total proton column density yields an elemental abundance. 
We adopted a combination of the H2 data of Savage et al. (1977) and H I data of Diplas & Savage 
(1994) for p Oph A, Nh = 5.0 x 10^^ cm-^, and the results of Bohlin, Savage, & Drake (1978) for 
X Oph, 2.26 X 10^^ cm~^. The fractional abundances (pa 10~^°) are similar to those seen toward 
C Oph (Cardelli, Savage, Sz Ebbets 1991; Cardelli ct al. 1993b) and in other H2-rich sight lines 
(Cardelh 1994 for Ga II; Sofia, Meyer, & Cardelh 1999 for Sn II). 

The archive of Copernicus data developed by George Sonneborn and available at the 
Multiwavelength Archive at the Space Telescope Science Institute was used in this research. 
Support for this work was provided by NASA through grant number GO-5389.02-93A from 
the Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Incorporated, under NASA contract NAS5-26555 and through Long Term 
Space Astrophysics grant NAG5-4957 to the University of Toledo. We acknowledge the helpful 
suggestions of an anonymous referee. 
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Table 1 
GHRS Observations 



Star 


Filename 


Wavelength 

(A) 


Counts 


Total Counts 


S/N Ratio 


CO bands 


p Oph A« 


z2b95106 


1465.60 


10547 










zzbyolUT 


140/ .04 


10363 










z2b95108 


1461.30 


10045 










z2b9510a 


1471.35 


6445 


37400 


193 


2-0, 3-0 




z2b9510b 


1403.67 


6983 










z2b9510c 


1405.12 


6489 










z2b9510e 


1408.00 


5989 


19461 


140 


4-0, 5-0 




z2b9510f 


1271.36 


8917 










z2b9510g 


1272.83 


8030 










z2b9510i 


1275.75 


7786 


24733 


157 


10-0, 11-0 




z2b90206 


1465.60 


15141 










z2b90207 


1467.04 


15388 










z2b90208 


1461.30 


15358 










z2b9020a 


1471.35 


14415 


60302 


246 


2-0, 3-0 




z2b9020b 


1403.67 


15559 










z2b9020c 


1405.12 


14975 










z2b9020e 


1408.00 


15279 


45813 


214 


4-0, 5-0 




z2b9020f 


1271.36 


22992 










z2b9020g 


1272.83 


17850 










z2b9020h 


1275.75 


23780 


64622 


254 


10-0, 11-0 



" GHRS exposures of HD147933 with G160M conducted on 1995 June 17. 
GHRS exposures of HD148184 with G160M conducted on 1995 February 11-12. 
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Table 2 
CO Measurements 



Molecule 



CO Band 



(mA) 









p Oph A 


X Oph 




A 


-X 2-0 


87.7 ±1.6 


48.9 ± 1.2 






3-0 


80.9 ± 2.0 


44.6 ± 1.2 






4-0 


62.5 ±2.0 


39.9 ± 1.2 






5-0 


53.3 ± 1.9 


28.3 ± 1.6 






10-0 


10.0 ± 1.4 


2.0 ±0.6 






11-0 


4.6 ± 1.6 


1.6 ±1.0 




o! - 


- X 14-0 


13.8 ± 1.5 


3.6 ±0.5 




e 


-X 5-0 


3.6 ±0.9 


0.8 ±0.6 




A 


-X 2-0 


10.3 ± 1.5 


2.6 ± 1.0 






3-0 


11.0 ± 1.6 


2.2 ± 1.1 






4-0 


4.9 ± 1.3 


1.5 ±0.8 






5-0 


4.1 ± 1.0 


0.8 ±1.0 
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Table 3 



Results for Vibrationally Excited H2 



Line 




Wx 
(mA) 






p Oph A " 


X Oph 


C Oph 


R(0) 


1.07 ±0.19 




0.34 ± 0.16 


R(l) 


2.09 ± 0.27 




0.33 ± 0.16 


R(2) 


2.00 ±0.35 


< 0.52 


< 0.30 


R(3) 


2.12 ±0.35 


< 0.68 




P(l) 


< 0.70 






P(2) 


0.81 ± 0.26 






P(3) 


1.43 ± 0.28 


< 0.53 





" Weighted average of independent determinations 
by J.A.C. and S.R.F. 

^ From Federman et al. (1995). 
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Table 4 





Results 


for and 


13QH+ a 




Star 


Molecule 


VFa 


VLSR 


6 






(mA) 


(km s^""^) 


(km s~^) 


p Oph A 




12.25 ±0.09 


4.5 


2.2 




13CH+ 


0.11 ±0.05 




2.2 




12CH+ 


10.45 ±0.16 


2.6 


2.1 






0.17 ±0.05 




2.1 



" Results based on separation of —0.272 A and the line width for ^^CH 
was constrained to be the same as that for ^^CH"*". 
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Table 5 

Results from Fitting the CO Bands 



Parameter p Opli A x Oph 
N(^'^CO) (em-2) (1.92 ±0.25) xlO^^ (3.8 ± 1.0) xlO^"^ 

?;/,eHo (km s-i) -4.4 ±1.6 -10.3 ±1.5 

b (km s-i) 0.60 ± 0.02 0.66 ± 0.04 

Tio (K) 2.7 ±0.1" 3.0 ±0.3^ 

T21 (K) 7.6 ± 0.5'^ 5.2 ± 0.6* 

T32 (K) 8.4 ± 0.5'^ 7.5 ± 3.7* 

iV(i2C0)/iV(i3C0) 125 ±23 117 ±35 

N{C}-^0)/N{C^^O) 1100 ±600 



" Tio, T21, T32 = 10.3 ± 2.9 K from the 10-0 and 11-0 bands, 
while the 2-0 to 5-0 bands of ^^CO give 7.7 ± 2.5 K. 
^10, T2U T32 = 4.7 ± 2.2 K from the 10-0 and 11-0 bands, 
while the 2-0 to 4-0 bands of ^^CO give 4.4 ± 2.9 K. 
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Table 6 

Comparison of Results for the B — X and C — X Bands 



Star 


Wb- 


-X (mA) 


Wc- 


-X (mA) 




Present 


Fit Other 


Present 


Fit Other 


p Oph A 


27.3 ±4.8 


32.4 


54.6 ±5.9 


65.1 22 ±7" 


X Oph 


9.6 ±3.7 


11.9'' 


42.1 ±4.1 


43.9 44 ± 2.5 " 



« Federman et aL 1980. 

Since the P branch is within the noise, the W\ for the R branch is given. 
" Frisch 1980. 
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Table 7 



Modeling Results for Vibrationally Excited H2 



Species 






Model 






Observed 


1 


2 


3 


4 


5 


p Oph A " 














luv 


4.0 


10.0 


10.0 


10.0 


8.0 




ntot (cm~^) 


300.0 


300.0 


300.0 


300.0'' 


300.0'' 




T(K) 


45.0 


100.0 


100.0 


100.0'' 


45.0'' 




Index 


-1.00 


-1.00 


-1.00 


-1.20 


-1.18 




Form, model 




ST 




ST 


ST 




H 


4.8(20)^ 


6.6(20) 


6.6(20) 


1.0(21) 


1.1(21) 


6.5(21) 


H2 


3.7(20) 


3.7(20) 


3.7(20) 


3.7(20) 


3.7(20) 


3.7(20) 


H2 u = 0: 














J = 


3.1(20) 


1.4(20) 


1.4(20) 


1.2(20) 


2.1(20) 


(3.0± 1.0)(20) 


J= 1 


6.4(19) 


2.3(20) 


2.3(20) 


2.5(20) 


1.6(20) 


(7.1 ±3.0)(19) 


J = 2 


1.7(17) 


2.7(18) 


2.7(18) 


4.9(18) 


1.6(18) 




J = 3 


1.1(16) 


5.0(16) 


5.0(16) 


1.7(17) 


7.9(16) 




J = 4 


1.5(15) 


2.1(15) 


2.3(15) 


2.3(15) 


1.8(15) 




J = 5 


2.9(14) 


7.0(14) 


7.2(14) 


7.7(14) 


5.2(14) 




H2 V = 3: 














J = 


5.4(11) 


7.0(11) 


7.0(11) 


6.8(11) 


5.6(11) 


(8.2±2.0)(11) 


J= 1 


7.4(11) 


2.3(12) 


2.3(12) 


2.6(12) 


1.6(12) 


(2.0±0.5)(12) 


J = 2 


1.3(12) 


1.7(12) 


1.7(12) 


1.7(12) 


1.4(12) 


(1.7±0.4)(12) 


J = 3 


6.4(11) 


2.0(12) 


1.9(12) 


2.2(12) 


1.4(12) 


(2.7±0.5)(12) 
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Table 7 Cont. 
Modeling Results for Vibrationally Excited H2 



Species 




Model 




Observed 


1/ 


2 a 


3^^ 


X Oph 










H2 u = 3: 










J = 


1.0(12) 


1.6(11) 


3.0(11) 




J = 1 


1.4(12) 


2.2(11) 


4.2(11) 




J = 2 


2.4(12) 


3.8(11) 


7.1(11) 


< 4.8(11) 


J = 3 


1.2(12) 


1.9(11) 


3.5(11) 


< 8.9(11) 



" All models have the standard extinction curve with grain model 2. 
^ Temperature and density gradient: values refer to center of cloud. 

ST: statistical distribution of 1.5 eV over all levels upon formation. 

W: Wagenblast 1992: fraction 0.73 in i; = 5, J = 9 and 0.27 in J = 10. 
^ 4.8(20) = 4.8 X 10^° cm-2. 

^ X Oph model G: T = 45 K, luv = 9, ritot = 300 cm'^ (van Dishoeck & Black 1986). 
3 Model G with luv = 1- 
^ Model G with luv = 2. 
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Table 8 



Results for Ga II, As II, and Sn II 



Species 


Line 


/-value 


p Oph A 


X Oph 




(A) 




Wx 


iV(X) (xlO^i) 


Wx 


N{X) (xlO^i) 








(mA) 


(cm"^) 


(niA) 


(cm~^) 


Ga II 


1414.40 


1.80 « 


11.28 ±0.35 


4.38 ±0.17 


0.81 ±0.16 


0.26 ± 0.05 ^ 












5.41 ±0.16 


1.87 ±0.06 ^ 


As II 


1263.77 


0.32 " 


1.23 ±0.35 


2.78 ±0.81 


0.67 ±0.24 


1.50 ±0.55 


Sn II 


1400.44 


1.04 


4.69 ± 0.33 


2.82 ± 0.22 


0.41 ±0.16 


0.24 ±0.09 












3.02 ±0.16 


1.76 ±0.10 ^ 



" Prom Morton 1991. 

There arc two velocity components, separated by approximately 10 km s 
" Prom Cardelh et al. 1993b. 
^ Prom Schectman et al. 2000. 
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Fig. 1. — a) Spectra of ^^CO absorption toward p Oph A. The data are represented by the filled 
circles. Our best fit to the data (solid line) and the data— fit (line, offset to 1.10) are also shown. 
The fit to the a' — X (14—0) intersystem band was not used in the derivation of the parameters 
given in Table 5. Absorption from As II is seen in the spectrum for the 11—0 band, (b) Spectra 
for X Oph. 
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Fig. 2. — a) Spectra of ^^CO absorption toward p Oph A. The data are represented by the fihed 
circles. Our best fit to the data (soHd hne) and the data— fit (hne, offset to 1.04) arc also shown. 
The fit to the e — X (5—0) intcrsystem band was not used in the derivation of the parameters given 
in Table 5. (b) Spectra for x Oph, except here the quality of the fit is offset to 1.03. 
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Fig. 2.— b) 
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Fig. 3. — Spectra showing absorption from vibrationally excited H2 toward p Oph A. The strong 
hnes are due to C I absorption. 
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Fig. 4. — Rectified spectra of CH+ absorption toward p Oph A and % Oph. The inserts focus on 
the region of ^^CH+ absorption. 
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Fig. 5. — Curves of growth for the A — X system of bands in ^^CO (solid Hncs) and ^^CO (dashed 
lines) toward p Oph A (filled circles) and x Oph. (open circles). The numbers along the top give v' 
for the band. The resulting ^^CO/^^CO ratios are also indicated. 
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Fig. 6. — Results of excitation analyses for C I and C2- The density indicated here refers to collision 
partners, ric- The dashed and dot-dashed curves show the constraints inferred from N{1)/N{0) 
and N{2)/N{0), respectively, while the solid contour represents the acceptable values from the C2 
rotational lines. The range in kinetic temperature deduced from the relative amounts of H2 in J = 
and 1 is indicated by dotted lines. For p Oph A, the gray curves show the effect of a larger UV 
flux on C I excitation. 



